Unification Theory of Angular Magnetoresistance Oscillations in 
Quasi-One-Dimensional Conductors 
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We present a unification theory of angular magnetoresistance oscillations, experimentally observed 
in quasi-one-dimensional organic conductors, by solving the Boltzmann kinetic equation in the ex- 
tended Brillouin zone. We find that, at commensurate directions of a magnetic field, resistivity 
exhibits strong minima. In two limiting cases, our general solution reduces to the results, previ- 
ously obtained for the Lebed Magic Angles and Lee-Naughton-Lebed oscillations. We demonstrate 
that our theoretical results are in good qualitative and quantitative agreement with the existing 
measurements of resistivity in (TMTSF)2C104 conductor. 

PACS numbers: 74.70.Kn, 72.15.Gd, 71.18.-|-y 



I. INTRODUCTION 

Magnetic properties of quasi-one-dimensional (QID) 
organic conductors have been intensively studied both ex- 
perimentally and theoretically since a discovery of the so- 
called Field-Induced Spin-Density- Wave (FISDW) phase 
diagrams^'^. For open electron orbits, where the Lan- 
dau levels quantization is impossible, as theoretically 
shown^''^"^, other quantum effects - the Bragg reflections 
of electrons from the Brillouin zones - play an important 
role. In the simplest situation, where magnetic field is 
perpendicular to conducting layers of QID conductors, 
the Bragg reflections are demonstrated^"^ to be responsi- 
ble for the FISDW phases formation. In an inclined mag- 
netic field, perpendicular to conducting chains, a more 
complicated interference picture of electrons, moving in 
the extended Brillouin zones, appears. As shown in Refs. 
6 and 7, it results in the constructive interference of elec- 
tron waves in some many-body effects for some special 
commensurate directions of a magnetic field, which are 
called the Lebed Magic Angles (LMA). 

The LMA effects were experimentally discovered in 
Refs. 8-13 and are observed in a number of organic 
conductors^^"^^, which possess QID parts in their Fermi 
surfaces (FS). Note that, instead of maxima of resistiv- 
ity due to electron-electron scattering, predicted in Ref. 
7, the experiments^"^^ demonstrate clear minima at the 
LMA directions of a magnetic field. In important theo- 
retical contributions^^'"^^, it was shown that constructive 
interference effects^'^ could appear in such one-body phe- 
nomenon as a residual resistivity due to impurities. As a 
result, minima of resistivity component, perpendicular to 
conducting layers, were theoretically found at the LMA 
directions of a magnetic field^°'^^. Nevertheless, theo- 
retical model^° predicted weak (i.e., exponentially small) 
magnitudes of the LMA minima, whereas Ref. 31 was 
based on an incorrect solution of the Boltzmann kinetic 
equation. 

Recently, a correct solution of the Boltzmann equation 
for a magnetic field, perpendicular to conducting chains 
of a QID conductor, was found and the existence of the 



strong LMA minima in perpendicular to conducting lay- 
ers component of resistivity was firmly established'^^. In 
addition, a quantum mechanical variant of the theory was 
suggested'^^ . Theory^'^ reveals quantum interference na- 
ture of periodic solutions of the Boltzmann kinetic equa- 
tion in the extended Brillouin zone"^^, which lead to the 
appearance of the LMA effects. According to Ref. 33, 
the LMA minima of resistivity appear due to changes of 
electron wave functions dimensionality from ID into 2D 
at commensurate directions of a magnetic field^'^ due 
to constructive interference effects. These interference 
effects appear between velocity component, perpendic- 
ular to conducting layers, and the density of states^'^. 
Note that the theory^^'^'^, which is a limiting case of the 
unification theory, suggested in this paper, can explain 
the experimental observations of the LMA effects only 



in resistivity 
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As to the observations of anomalously 



strong LMA phenomena in the Nernst'^*"'^* and HalP^ 
effects, their explanations may need a different theoreti- 
cal approach. 

As experimentally discovered''""^'^, the LMA- like mag- 
netoresistance minima with enhanced magnitudes are 
observed^'^"^® in experimental geometry, where magnetic 
field is inclined with respect to conducting chains of 
a QID conductor. At first, these angular oscillations 
were interpreted^^'''^^ in terms of the LMA phenomenon. 
Later, it became clear'^^"^" that, although they are re- 
lated to the LMA effects, their concrete physical meaning 
is quite different. Below, we call the above mentioned an- 
gular oscillations of resistivity the Lee-Naughton-Lebed 
(LNL) ones. 

Theory of the LNL phenomenon, based on a solution of 
the Boltzmann kinetic equation in the extended Brillouin 
zone, was suggested in Refs. 47, 48. In Refs. 47, 49, the- 
ory of the LNL oscillations was extended to the so-called 
weak non Fermi-liquid quantum case^^. In Refs. 44, 
45, 50, the quantum theory of the LNL oscillations was 
suggested for the Fermi-liquid case and their quantum 
interference nature was revealed. According to Refs. 44, 
45, 50, the LNL angular oscillations are due to changes of 
electron wave functions dimensionality from ID into 2D 
at some commensurate directions of a magnetic field due 
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to constructive interference effects. In contrast to case 
of the LMA oscillations, these interference effects appear 
between two velocity components, perpendicular to the 
conducting chains. By present time, the physical origin of 
the LNL oscillations has been firmly established"'^~'^° and 
a comparison of the theoretical results with the existing 
experimental data have been made^'''^^'^''. These include 
theories of the LNL phenomenon in the presence of an- 
ion ordering potentials and theory^^, connecting 
the LNL oscillations with the Aharonov-Bohm effect. 

The goal of our paper is to suggest an analytical unified 
theory, which describes both the LMA and LNL phenom- 
ena as its limiting cases. Note that the suggested theory 
also describes the so-called Danner-Kang-Chaikin (DKC) 
osciUations^^ and Third Angular Effect (TAE)54-56, To 
derive an analytical expression for resistivity component, 
perpendicular to conducting layers, we analytically solve 
the Boltzmann kinetic equation in the extended Brillouin 
zone in a magnetic field, inclined with respect to conduct- 
ing chains of a QID conductor. We find strong minima 
of resistivity, corresponding to the appearance of peri- 
odic solutions of the Boltzmann equation at the LMA 
and LNL commensurate directions of a magnetic field. 
A comparison of our theoretical results with the experi- 
mental data on (TMTSF)2C104 conductor^^ shows good 
qualitative and quantitative agreement (see Figs. [Tl[2|). 



II. UNIFICATION THEORY 

Let us consider a QID conductor with the following 
electron spectrum, 

e=^(p) = ±Vrc{pyb*)[px :^p^[py)] - 2tcC0s(pzC*), 
Px{Py) ^Pf - 2tbg{pyb*)/vF. (1) 

[Here, +(— ) stands for right (left) sheet of the FS; 
vp and pf are the Fermi velocity and Fermi momen- 
tum along the most conducting x axis, respectively; tb 
and tc are the hopping integrals along y and z axes.] 
Note that in Eq. ([T|), in accordance with Refs. 31-33, 
Py dependence of the velocity along conducting chains, 
Vx{Py), is taken into account. For most QID conduc- 
tors, we can choose g{pyb*) — cos{pyb*) in Eq. ([!]), but 
for an important exception of (TMTSF)2C104 conduc- 
tor with an anion ordering, we have to use g{pyb*) = 
[cos^{pyb*) + (A/2t,,)2]i/2^ where A is the the so-called 

1 2 44 4^ 

anion gap ' ' ' . 

Below, we show that, when a QID conductor (IT|) is 
placed in a tilted magnetic field. 



H = iJ (cos 9 cos (f, cos 6* sin (p, sin 9) , 



(2) 



and a weak electric field perpendicular to the conducting 
(x, y) plane, then at certain orientations of the field. 



(3) 



sm tp — n [ — ] tan ( 



where n is an integer, resistivity, pzz{H,9,(p), exhibits 
strong minima. 

The Boltzmann kinetic equation in crossed electric and 
magnetic fields in the so-called r-approximation can be 
written in a standard way. 



eE 



[v X H] 



dfjp) _ /(p)-/o(p) 
dp T 



(4) 



[In Eq. (U), /(p) is an electron distribution function and 
T is a relaxation time.] After the standard approxima- 
tion, 

/(P) = /0(P) - ^^^{Py,P.), (5) 

the Boltzmann equation ^ can be written as^®, 

eE.v~f^)[vxH].^%i^ = ^^(^, (6) 
Vc/ op T 

where /o(p) is the Fermi-Dirac distribution function. 
[Note that, in Eq.(6), we use independent variables 
{e,Py,Pz)^^ , instead of {px,Py,Pz), where energy e is con- 
served in the absence of electric field. Since both electric 
field and temperature are supposed to be small, the elec- 
tric conductivity is defined by electrons, located in the 
near vicinity of the QID FS. Therefore, the distribution 
function '^{py,pz) can be taken at e = Ci? and, thus, 
does not depend on energy in Eq.(6).] Taking into ac- 
count that V — i9e(p)/c}p, we can rewrite the Boltzmann 
equation as, 

eEv^^ sin z - .,(,, 9)'-^ - ^, ^) 



dy 

dz T 



(7) 



To simplify the notations, we define the following dimen- 
sionless parameters: y = Pyb* , z = PzC* , and the follow- 
ing frequency variables: 



My,0) = [^yAy)Hb*sm9, 
Lu*{0,(p) — (^-^VyHc* COS 9 COS ip, 



cos f sm ip, 



(8) 



where = 2tbb* and v° = 2tcC*; h^l. 

It is important that the partial differential equation ([7]) 
can be analytically solved (see the Appendix A), 



.0 

<f{y,z) = eEv° / 

J —c 



du 



-oo ujb{y + u,0) 



dv 



In ^^biy + v, 
+<{o,v)g' {y + v) 

dv 



■ exp 



„ TUb{y + v,t 



(9) 
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Since the solution ([9]) of the Boltzmann equation ([7]) is 
known, we can express electric current density along z 
axis, perpendicular to conducting layers, in the following 

way, 



Pzz(0)(arb. units) 



dy 



dzz;0sinz^4^vl/(y,z) 
oe 



(10) 



In Eq. (jlOl) . we omit an exact factor since conductivity 
in a magnetic field is scaled below to its value in the 
absence of the field. Notice that dfo{p)/de in Eq. PU)) 
is the density of states at the FS, which is proportional 
to l/vx{y). From Eqs. ©, (ITU|) . we obtain the following 
result for inter-layer conductivity after some calculations. 



■ cos 



„ ujb{y + v,e) 



dy r du 

vx{y) 7-00 •^b{y + u) 

ujc{y + v,9,(p) 



■ exp 



d?; 



ujbiy + v)i 



(11) 



[Note that, in a QID case resistivity component, per- 
pendicular to conducting layers, is pzz{H,9,ip) — 
l/azz{H,9,ip).] Eqs.(9)-(11) and comparison of Eq.(ll) 
limiting cases with the experimental data, obtained on 
(TMTSF)2C104 conductor^''', are the main results of our 
paper. 



III. LIMITING CASES 

In this section, we consider two different experimental 
settings, which correspond to the LNL and LMA phe- 
nomena, repectively. 



A. LNL Oscillations 

In experiments, where the LNL oscillations are ob- 
served, magnetic field is tilted away from (y, z) plane. In 
the case, where 9 <^ 90°, the py dependence of the Fermi 
velocity is not significant, and, therefore, we can every- 
where replace Vx{py) into vp- This significantly simplifies 
the formula for conductivity dill) . 

dzexp(z) / — cos {uJc{9,(p)tz 
-oo Jo ^'^ ^ 



0Jb{9) 



(12) 



where the frequency variables are defined in Eq. ([5]) with 
''^xiy) = vp. The existence of resistivity minima at the 
commensurate directions of a magnetic field Q can be 
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FIG. 1: Comparison of the theory p2[) and the 
experimenl)^ for the Lee-Naughton-Lebed (LNL) oscillations 
in (TMTSF)2C104. Thin solid curve is calculated from Eq. 
(|12|) aX 9 = 7° , with band parameters being ta/tb = 8.5 and 
A/2tb — 0.1. For calculations, we have used ujb{0)T = 15 and 
c* = 26* . Numerical results are in good qualitative and quan- 
titative agreement with the experimental data (thick solid 
line) in a broad range of angle (f>. 



understood from Eq. (fT^ . Since g(jj) is a periodic func- 
tion of y, at the commensurate directions [where a con- 
dition ujc{9,ip) = nuJb{9) is satisfied], constructive inter- 
ference effects increase the integral p2|) . which gives the 
minima in pzz{H, 9, Lp). 

In Fig. [1] we show both the theoretical re- 
sults, obtained by means of numerical calculations of 
Eq. (fT2|) . and the experimental data^*", obtained on 
(TMTSF)2C104 conductor. Note that the minima of 
resistivity appear only for even values of n in Eq. 
([3]) due to the existence of the anion ordering gap, 
A, in (TMTSF)2C104 electron spectrum, g{pyb*) = 
[cos'' {pyh*) + (A/24)2]i/2 [see Eq.©]. 

Most QID conductors [e.g., (TMTSF)2PF6 and k- 
(ET)2Cu(NCS)2] can be described by the electronic spec- 
trum ID) with g{y) = cosy. In this case, Eq. reduces 
to, 

/■° dy r 

azz{H,9,Lp) ^ I d2exp(z) / —cos uJc{9,(p)tz 



M9) 
Using the identity. 



2tt 



(cos(?/) - cos (y + uJbi9)Tz)) 



+ 00 



(13) 



exp(i2Cos6') = ^ Jn{z)i'^e^ 



(14) 



n— — oo 



we can further simplify Eq. (|13p. 



<^zziH,9,(f) ^ dzexp(2:) 

^ — OO 



(15) 



+ 00 



ujt{9) 



cos {{uJc{9, (f) - nujb{9))Tz) 
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Finally, Eq. p5|) can be transformed into a simple ana- 
lytical form, 



Pzz{ci')(arb. units) 



a..(i/ = 0) 



n— — oo 



where Jn{---) is the n-oder Bessel function. The oscil- 
latory behavior of interlayer resistivity is directly seen 
from Eq. ()16p . Indeed, at the commensurate directions 
of a magnetic field ([3]), where ujc{0,(p) — nuji,{d), there 
appear maxima of conductivity (jl6p . which lead to re- 
sistivity minima in pzz{H,9,ip) (see Fig. [ij. Note that 
Eq.(16) is an agreement with the previous results [47-50]. 



B. LMA Effects 

The LMA phenomena are experimentally observed in 
a magnetic field, directed in (y, z) plane. 



H — {0, H sin a, H cos a) , 



(17) 



where a = 90° — 9. Therefore, Eq. ^ for the commen- 
surate directions of a magnetic field reduces to. 



tan a ~ n 



(18) 



where n is integer. In this case, where (f> = 90", Eq. (fTTj) 
for interlayer conductivity can be rewritten as. 



dy 


■f 




J-oo 


exp 


i-f 



du 



cos [N{a)u] 
ujbiy + u,a) 
dv 



„ TU}b{y + v,a) 



, (19) 



where ujb{y,a) = eHvx{y)b* cos a/ c, tudy^a) — 
eHvxiy)c* sina/c, and N{a) = iOc{y,a)/ujb{y,a). Eq. 
(|19p can be transformed to the following expression by 
integrating by parts. 



azz{H,a) j - 

J V: 



dy 



(y) 

• exp 



1 + N{a) / du sin [N{a)u 
dv 



„ TUJb{y + v,a) 



(20) 



If we introduce the following notations, 

VF 



f{y) 



vx{y) 

e 



1, 



hb{H,a) — -HvFb*T cosa, 
c 

hc{H,a) — -Hvpc*T sina, 
c 



(21) 



where l/vp = (l/uj.(y))j,^, then Eq. (1^ becomes, 
a.AO)- I dy{l + J{y)) 




Q'(degree) 



FIG. 2: Comparison of the theory (I26|l and experiment^ 
for the Lebed Magic Angle (LMA) phenomena in 
(TMTSF)2C104, where g(pyb') = [cos^{j)yh*) + {l^/2tbfY^^ 
in Eq.([l}. Thin solid curve is numerically calculated from Eq. 
(|26p . with band parameters ta/th = 8.5 and A/2ti, = 0.1. For 
the calculations, we have used cji,(0)r = 15 and c* = 2b*. 
There exist good qualitative and quantitative agreement be- 
tween the theory and experiment for 0° < a < 60°. For 
negative values of a, agreement between the theory and ex- 
periment is worse, perhaps, due to some problems with the 
experimental measurements. 



l + hc{H,a) I du sin [hc{H, 

J —oo 



a)u] 



■ exp 



dvil + f{y + hbiH,a)v)) 



(22) 



Integrating by parts one more time and taking into ac- 
count a periodicity of function f{y), we obtain the fol- 
lowing expression for conductivity. 



"'f^'^l ^l-hliH^a) f d.exp 

(Tzz[tl = 0) 



(u) 



(23) 



/■^'^ dy r f° 
■cos[hc{H,a)u] J 27r^^^| ./ ^"^^fiy + ^^iH, 



a)ui 



In the so-called clean limit, where hb{H,a) 1, Eq. 
([23l) can be significantly simplified. Below, we introduce 
the Fourier transform of function f{y), 



+ 00 



f{y) = ^Ancosiny). 



(24) 



In the clean limit, the last exponential function in Eq. 
(1^51) . whose argument is inversely proportional to hb{H), 
can be expanded as, 

expl^-J duif[y + hb{H,a)ui]^ = 1 

dui ^ An cos [n{y + hb{H, a)ui)] 

-f- / dui / du2 ^ An Am COS [n{y + hb{H,a)ui)] 



■ COS [m{y + hb{H, a)u2)] , 



(25) 
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where the higher order terms are discarded. After in- 
tegration with respect to variable y, the second term in 
Eq. ()25|) vanishes, whereas, in the third term, only con- 
tributions with n — m retain. Finally, the interlayer 
conductivity can be represented as, 



azziH, a) 



1 



a,,{H = 0) l + hl{H,a) 

^Al( 2 



tan a 



2b* 



1 



1 -I- [hc{H,a) - nhb{H,a)Y 
1 

l + [hc{H,a)+nhb{H,a)Y 



(26) 



Note that Eq.(26) is an agreement with our previous re- 
sults [32,33]. 

In Fig. [21 we compare the theory (pS)) with the exper- 
imental data^^ for the LMA phenomenon. It is impor- 
tant that, for the calculations, we have used in Eq. ipS)) 
the same values of the parameters as for the calculations 
of the LNL phenomenon [see Eq. (fT2|) and Fig. [T]. Fig. 
[2] demonstrates good qualitative and quantitative agree- 
ment between the theory and experiment in the wide 
range of the angles: 0° < a < 60°. For |a| > 60°, 
there appear significant deviations from the experimental 
behavior^^. One possible reason for that is a breakdown 
of Ea. (p6)) for large values of angle a, where the so-called 
clean limit approximation is not valid. Another possible 
reason for the deviations is that, at high values of angle 
a (i.e., at high in-plane projections of a magnetic field), 
there may occur Fcrmi-liquid'^^ or non Fcrmi-liquid^*^ de- 
coupling of the conducting layers, where the Boltzmann 
kinetic equation is not valid any more. 



IV. CONCLUSION 

In this work, we propose a unification theory for an- 
gular magnetoresistance oscillations, experimentally ob- 
served in QID organic conductors. We analytically solve 
the Boltzmann equation in the extended Brillouin zone 
and find analytical formula, which describes interlayer re- 
sistivity. We show that, in two important limiting cases, 
this formula reduces to the expressions, previously ob- 
tained to describe the LNL and LMA phenomena in re- 
sistivity. Numerical results, obtained from these expres- 
sions, are shown to be in good agreement with the exper- 
imental data, obtained on (TMTSF)2C104 conductor, in 
a broad rage of magnetic field directions. On the other 
hand, a comparison of the theory with the LMA exper- 
imental data reveals significant discrepancy between the 
theory and experiment for directions of a magnetic field 
close to the conducting layers. We suggest that this dis- 
crepancy may be due to decoupling of the conducting 
layers in a parallel magnetic field. 
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Appendix A: Proof of Eq. ((9|) 

From physical point of view, we have to find a par- 
ticular solution of the inhomogeneous partial differential 
equation (7), which is proportional to external electric 
field and periodic in the extended Brillouin zone with re- 
spect to both variables, y and z. To find such solution, 
below we use the method of characteristics for the first 
order partial differential Eq. (7). Let us consider the 
following equation, which defines the characteristics: 



dz ^ Wcjy , 0, (f) + u*{d, (p)g'{y) 
dy 



(Al) 



^b{y,0) 

with initial condition being, 

z{u) = ZQ. (A2) 
Note that Eqs. (Al),(A2) have the following solution, 
'■^ dv 



z{y) = zq + 



[Lu,{v,e,'f)+ujl{e,'f)g'{v)] 



.(A3) 

At this point, Eq. ^ can be considered as an ordinary 
differential equation, 

eEv'!. sm zo ~ujb[y, 0) = . ( A4) 

dy T 

For our purpose, we choose the following solution of Eq. 
(A4), 



J — ( 



dv 



ujb{v,6) 



sm zq exp 



dw 



/„ UJb{w,9)T 

(A5) 

If we take into account that zo is given by Eq.(A3), we 
can obtain: 



■^{y,z) = eEvl T 

J — c 



du 



-oo ^b{u 

^ dv 



u Wfc(w, 



+u^l{6,^)g' (v) 



exp 



dv 



TUJb{v,0) 



(A6) 



From Eq. (A6), it is directly seen that the obtained solu- 
tion of Eq.(7) is proportional to electric field and periodic 
with respect to variable z. It is possible to make sure 
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that it is also periodic with respect to variable y. For 
this purpose, it is necessary to take into account that 
the functions uJb{y,0), uJc{y,0,(j)), and g{y) are periodic 
with respect to variable y in the extended Brillouin zone. 



After changes of integration variables, u ^ y + u and 
V ^ y + V, we can finally obtain Eq. (9) from Eq. (A6). 

* Also at: Landau Institute for Theoretical Physics, 2 
Kosygina Street, 117334 Moscow, Russia. 
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